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ABSTRACT 



Aims. Combining the relative vicinity of the Local Group spiral galaxy M 33 with the Spitzer images, we investigate the properties of 
infrared (IR) emission sites and assess the reliability of the IR emission as a star formation tracer 

Methods. The mid- and far-IR emission of M 33 was obtained from IRAC and MIPS images from the Spitzer archive. We compared 
the photometric results for several samples of three known types of discrete sources (Hn regions, supernovae remnants and planetary 
nebulae) with theoretical diagnostic diagrams, and derived the spectral energy distribution (from 3.6 to 24 fim) of each type of object. 
Moreover, we generated a catalogue of 24 pm sources and inferred their nature from the observed and theoretical colours of the known 
type sources. We estimated the star formation rate in M 33 both globally and locally, from the IR emission and from the Ha emission 
I ■ line. 

Results. The colours of the typical IR emissions of Hn regions, supernovae remnants and planetary nebulae are continuous among 
the different samples, with overlapping regions in the diagnostic diagrams. The comparison between the model results and the colours 
of Hn regions indicates a dusty envelope at relatively high temperatures ~600 K, and moderate extinction Ay^ 10. The 24 pm sources 
IR colours follow the regions observationally defined by the three classes of known objects but the majority of them represent Hn 
regions. The derived total IR luminosity function is in fact very similar to the Hn luminosity function observed in the Milky Way and 
in other late type spirals. Even though our completeness limit is 5 x 10 37 ergs s _1 , in low density regions we are able to detect sources 
five times fainter than this, corresponding to the faintest possible Hn region. The 8 and 24 pm luminosities within the central 5 kpc of 
M33 are comparable and of order 4 x 10 28 ergs s" 1 Hz" 1 (vL v (8) = 1.5 x 10 42 and vL v (24) = 4.4 x 10 41 ergs s" 1 ). We estimate the 
total IR emission in the same region of M 33 to be 10 9 L Q . The discrete sources account for about one third of the 24 pm emission 
while the rest is diffuse. From the IR emission, we derive a star formation rate for the inner disk equal to 0.2 M© yr _1 , consistent with 
the star formation rate obtained from the Ha emission. 
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■ 1. Introduction In this series of papers, we examine star formation proper- 
ties of the Local Group late type spiral M 33. The proximity of 
Local Group members together with Spitzer resolution makes 

X " The interstellar medium (ISM) in galaxies provides the raw ma- possible the resolution of structures on spatial scales of 5-10 pc, 

^ terial from which stars form, and is the repository of the products and enable us to effectively study individual star forming sites 

■ - ' of stellar evolution. Stars form in condensations of cool molecu- and their surrounding environment in a galaxy different from our 

lar gas, eventually destroying their birth sites through energetic own. M 33, is at a distance of 840 kpc (|Freedman et al|1221l), its 

photons, massive stellar winds, and supernova (SN) explosions. mass and apparent size are smaller than that of M 3 1 but M 33 

Dust is formed through coagulation in stellar winds and super- hosts the brightest Hn complex in the Local Group. This, to- 

nova ejecta, and over time, becomes an integral part of the ISM. gather with its blue colour assures us that star formation is still 

Radiative and mechanical energy are input to the ISM through very active throughout its disk. M33 bears no signs of recent 

massive star formation and evolution, in the form of ionising mergers and therefore gives us a unique opportunity for investi- 

photons, supernova explosions, and magnetic fields. The ISM gating the interplay of gas, dust, and star formation in isolated 

properties, such as magnetic field, gas, and dust, affect star for- galaxy disks. This is possible thanks also to high resolution and 

mation processes on all scales. At the same time, star forma- sensitive observations now becoming available for this nearby 

tion drives the evolution of the ISM. On the long term, local galaxy at all wavelengths. Chemical abundances h ave been mea- 

processes of star formation can influence large-scale structure sured in stellar populations of different ages (e.g. |Magrini etal] 

in the ISM and determine how galaxies evolve. Ultraviolet and |2007afl), and surveys of atomic and m olecular hydrogen point 

optical wavelengths are not the best wavelengths for SF stud- out the location o f massive gas clumps ([En gar giola etajj [20031; 

ies, because dust can hide massive star formation, mask the |Heyer et alJlZQQfl). Through detailed dynamical mass modeling 

effects of feedback, and re distribute the ISM energ y and radi- of the galaxy we can now trace the mass distribution of visible 

ation. However, firs t IRAS dNeugebauCTetZIH98^. then ISO and dark matter (|Corbelli & Salucci||2000fc |Corbelli||20Q3). 

(Kess ler et al.lll996l) . and now Spitzer (iWerner et al.ll2004l) have Our focus here is first on M33's population of discrete 

been able to penetrate the dust in the ISM, and thus better study known type of point sources: Hn regions, supernova remnants 

the interaction among its constituents. (SNRs), and planetary nebulae (PNe). In principle, these sources 
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track different phases of the star formation process, from the 
generation of ionising photons and stellar winds in massive stars, 
to subsequent explosion as SNe and finally to the more evolved 
stage of PNe. Hence, our immediate goal is to assess the physi- 
cal conditions of the stellar environment at different stages. The 
second purpose of this paper is the generation of an IR source 
catalogue from the Spitzer MIPS observations at 24 pm. Here 
we shall use IRAC-MIPS diagnostic diagrams of known sources 
to better understand the nature of these sources. In a subsequent 
paper we will relate the IR properties with detailed observations 
of the surrounding ISM at other wavelengths and to the disk 
large-scale structure. Spitzer data and the vast multi- wavelength 
dataset available for M 33 will help understand how star forma- 
tion proceeds in low luminosity spiral galaxies and how this re- 
lates to other global galaxy properties. 

The present paper is structured as follows: in Sect. [2 we 
present the Spitzer IRAC and MIPS data and their reduction pro- 
cess. The large scale structure of the dust emission in M 33 is 
investigated through colour images in Sect. [3] We analyse the 
IRAC and MIPS 24 photometry of discrete type-known sources 
(Hn regions, PNe, SNRs) in Sect.@]and compare their colours to 
theoretical diagnostic diagrams in Sect. [5] In Sect. we present 
a catalogue of the sources emitting in the 24 pm MIPS band 
and interpret their nature in light of the results obtained in the 
previous Section. The reliability of the IR emission as a star for- 
mation tracer, through comparison with Ha, is investigated in 
Sect. [71 Finally, the summary and conclusions of our study are 
given in Sect. 



2. Observations and data reduction 

We retrieved images from the Guaranteed Time Observations 
(PIP 5, PI R. G ehrz) in the Spitzer Space Telescope 
(iWerner et al.ll2004l) data archive. The Spitzer Space Telescope 
carries two photometric cameras: the Infrared Array Camera 
(IRAC) and the Multiband Imaging Photometer (MIPS). 

2.1. IRAC data 

The Infrared Array Camera (iFazio et al.ll2004|) is equipped with 
two detector arrays of 256x256 pixels. The field of view of a sin- 
gle array is 5 .'2 x 5 .'2. We analyze here eight sets of IRAC obser- 
vations of M33 (AORs 3636224, 3636480, 3637760, 3638016, 
3638784, 3639040, 3640320, 3640576) in all four IRAC bands. 
The Basic Calibrated Data (BCD) were created by the Spitzer 
Science Center (SSC) pipeline, version S 14.0.0 for all AORs. As 
the first data frame of each observation sequence has a shorter 
integration time than the regular exposure time (10.4 s), we dis- 
carded those frames. A more complet e description of the ob - 
ser vations is given bylHinz et al.l (l2004l) . lMcOuinn et al.l (120071) . 
and lTabatabaei et al.l (l2007ah . 

The mosaics were assembled gathering all the BCDs for each 
wavelength. The individual calibrated frames were processed us- 
ing Mopex (iMakovoz & Khan 2005), with a cosmic-ray rejec- 
tion and a background matching applied between overlapping 
fields of view. During the reduction, we used a common value 
of the pixel size equal to 1720. The final number of individual 
BCDs used is about 1700, in each IRAC channel. The final di- 
mension of the mosaics is approximately 62.'4 x 91. '7. Because 
of the overlap between adjacent BCDs, the mean redundancy is 
12 (from 6 to 30). 

The obtained IRAC mosaics were background subtracted, 
estimating 20 sky levels (the median values in 5 x 5 pixels boxes) 



near the edges of the mosaics, the farthest from the galaxy cen- 
tre. The 5.8 and 8.0 pm channels presented a gradient in the 
background, and we removed it using a fitted surface with the 
IRAF task imsurfit (images. imfit). In each channel, the final 
background value subtracted is the mean of the 20 median val- 
ues: 0.059 at 3.6 pm, 0.146 at 4.5 pm, 1.475 at 5.8 pm, and 5.389 
at 8.0 pm. The mosaics were then aligned using the IRAF tasks 
geomap and geotran (in immatch) using 15 point sources. 

Because of the diffraction limit of the telescope, the images 
have a resolution of 177 to 270 (see IRAC Observing Manual) 
which translate into 6.7 to 8 pc at the distance of M 33. The final 
IRAC mosaics at 3.6, 4.5, 5.8, and 8.0 pm are shown in Fig.Q] 

2.2. MIPS data 

Images of M33 at wavelengths 24, 70, and 160 pm (with 
bandwidths of 5, 19, and 35 pm, respectively) were obtained 
using data collec ted by the Multiband Imaging Photometer 
dRiekeet al.ll2004l) . MIPS uses true detector arrays: 128 xl28 
pixels at 24 pm (Si:As detector), 32 x 32 pixels at 70 pm (Ge:Ga 
device) and 2 x 20 pixels at 160 pm (Ge:Ga device). The nom- 
inal fields of view are 5' x 5' at 24 and 70 pm and 0.'75 x 5' at 
160 pm. The resolution of the instrument is diffraction limited 
to 6", 18", and 40" at 24, 70, and 160 pm, respectively. The 
operating modes of all the observations were "scan map", the 
telescope scanning the sky at a constant rate. 

As for IRAC, we retrieved eight MIPS datasets from the 
archive (AORs 3647744, 3648000, 3648256, 3648512, 3648768, 
3649024, 3649280 and 3649536). The 24 pm BCDs were cre- 
ated by the SSC pipeline version SI 6. 1.0. The exposure times of 
the DCEs were 3.67 s (except for the first frames of each series, 
where the integration time per pixel is lower and which conse- 
quently were not used in the further reduction steps). Since the 
observations were acquired over two epochs, the respective mo- 
saics were assembled separately using Mopex. For consistence 
with the IRAC reduction, the pixel size was chosen to be 1720. 
The first epoch includes 4199 frames, and the second one in- 
cludes 4500 frames. This latter map is displayed in Fig. [2] The 
mean redundancy per sky position is 22 (varying from 10 to 40), 
in each mosaic. This high level of redundancy removes any spu- 
rious effects such as cosmic-ray hits and provides more reliable 
measurements. In the present article, as we plan to focus on point 
sources (Hn regions, SNRs and PNe) located near the centre of 
the galaxy, we chose to concentrate, for each epoch, on the two 
central AORs. These two images were then aligned with the ref- 
erence of the coordinates of 8 point sources and combined to- 
gether by keeping the average of the pixel values. A small back- 
ground gradient was still visible in the combined image and was 
removed by fitting a surface to selected regions (IRAF task: im- 
surfit). The zone defined by the superposition of the two central 
AORs of each epoch is marked on the 24 pm image (see Fig. [2]). 
The resolution measured on the final image is 6". 

The 70 and 160 pm MIPS observations were acquired with 
the same AORs as the 24 pm observations. The version SI 6. 1.0 
of the SSC pipeline was also used to create the BCDs. The inte- 
gration times of the DCEs is 4.19 s, both for the 70 and 160 pm 
observations. We used Mopex to create the final mosaics (again 
the pixel sizes were chosen to be 172), including 8750 individ- 
ual frames at each of the two wavelength channels. The mean 
redundancy is about 20 (from 8 to 42) and 6 (from 2 to 10) at 
70 and 160 pm, respectively. The final images show some strip- 
ing along the scan direction that is a residual instrumental arti- 
facts du^tojhetim^ of the Ge:Ga detec- 
tors dHelou et al.ll2004l: Iffinz et al.ll2004l) . The final resolution is 
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Fig. 2. MIPS images of M 33. The upper panels display the 24 fim images: the full image (one single epoch) is on the left. The inner 
region enclosed by the continuous line will be further investigated in the present paper, and is shown in the upper right panel after 
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about 16" and 40" for the 70 and 160 /mi images, respectively. 
The final MIPS images (24, 70, and 160 /mi) are shown in Fig. [2] 

2.3. Ha image 

We also include in our analysis an optical narrow-band Ha im- 
age, a standard star formation indicator. The M 33 Ha emission 
was observed at the Kitt Peak National Observatory (NOAO). 
The dimensions of the CCD are 2048 x 2048 with pixels of 
2! 03 and a total field of view of about 70' x 70'. More de- 
tails of the observations an d reduction process can be found in 
iHoopes & Wa lterbos ( 2000); the estimated [N n] contamination 
accounts for at most 5% of the measured flux. The total Ha lumi- 
nosity of M 33 is abou t 7 x 10 6 L o and is dominated by Hn regions 
(iDevereux et al.l 19971) . The optical extinction in M 33 is found to 
be globally low, A v < 0.4 mag. Towards Hn regions Ay~ 1 mag 
on average even though it can reach values as high as 2 mag 
towards high surface brigh tness cores (e.g. Isra el & Kennicuttl 
ll980tlDevereux et alii 19971) . The moderate values of extinction 
imply that Ha can be used as an SF tracer globally, but towards 
individual sources an extinction correction can become neces- 
sary. 

3. Large scale structure of stars and dust 

In the present paper, we concentrate on the inner region of M 33 
as defined by the Spitzer inner AORs which is outlined by the 
continuous line in Figs. [T] and [2] for the IRAC and MIPS data, 
respectively. 

In order to help the interpretation of Spitzer data, we first re- 
view briefly which are the main emission phenomena that dom- 
inate each IRAC and MIPS waveband. 

3.1. Origin of the emission detected by Spitzer 

The stars (mainly Red Giant Branch and/or Asymptotic Giant 
Branch, according to the age of the population) account for most 
of the emission in the 3.6 /mi IRAC band. A mix of stars and 
very hot dust continuum contribute to the emission in the 4.5 //m 
band. Nevertheless, in discrete sources with negligible contin- 
uum (see Sect. 0]), the emission seen in both these bands can 
have a strong component of ionic and molecular lines (e.g., 
iReach et al.l I2006I: IWilliams et aTll2006l) . In the 5.8 /mi IRAC 
band, the stellar contribution decreases and the dust continuum 
component gets stronger. The longest- wavelength 8.0 //m IRAC 
band still traces hot dust continuum but in most cases, the poly- 
cyclic aromatic hydrocarbons (PAH) dominates. If there are no 
PAHs and low continuum, however, in discrete sources there can 
be a sizable contribution fr om two H2 lines and the [Ar III] at 
8.9 //m (e.g. jHora etap2004l) . 

The PAH emission features (mainly concentrated at 3.3, 6.2, 
7.7, 8.6, and 11.3 iim ) are the optically active vibrational modes 
of PAH molecules dLeger & Pugetl 1 19841 lAUamandola et al.l 
fl989l) . resulting from internal conversions of energy following 
absorption of an optical or UV photon. The IRAC bands 3.6, 5.8, 
and 8.0 /mi contain P AH emission features at 3.3, 6.2, and 7.7 & 
8.6 /mi, respectively (lDraindl2003l:IChurchwell et al.ll2004l) . The 
4.5 jt/m channel is the only one with no PAH features, although it 
contains the bright hydrogen recombination line Bra (4.05 /mi). 

We investigate the far-infrared (FIR) emission with the MIPS 
images. The 24 /mi emission is contributed partly by thermal 
emission of classic (large) dust grains and partly b y very small 
grains (VSGs). Indeed, VSGs were introduced by iDesert et al.l 



(1990) as the carrier of the emission in the 25 /mi IRAS band. 
The 70 fim emission traces cooler dust, and the 160 //m one 
shows the diffuse cold dust component. 

3.2. Structure of infrared emission in M33 

Figure[T]shows that the stellar component is smoothly distributed 
everywhere in the disk of M 33. The warm and hot dust traced by 
the 5.8 /mi channel smoothly follow the morphological patterns 
(spiral arms and flocculent structures) of M 33 with a higher con- 
centration in the inner part of the galaxy. Compared with the 
other components, the 8.0 /mi contribution shows most clumpy 
distribution; evidently PAHs are preferentially located along the 
spiral arms and flocculent structures, rather than being uniformly 
distributed in the ISM. The very centre of the galaxy gathers all 
the three contributions: stars, dust and PAHs. 

3.3. The colour images 

To help investigate the large scale structure of M 33 as seen by 
Spitzer we made two sets of colour images: one degraded to 
the worst resolution of the IRAC channels, and a second set 
degraded to the resolution of the MIPS 24 /mi image. We ap- 
proximated the PRF with a Gaussian and convolved both sets to 
the resolution of the worst image in the set (8 /mi and 24 /mi, 
respectively). 

Fig. [3] shows the 3.6/4.5 ratio (left panel) and 4.5/8.0 ratio 
(right) of M33. The 3.6/4.5 shows a remarkably constant ratio 
(~ 1.7) over almost t he entire galaxy, confirming the statement 
by iPahre et al.l (|2004|) . who find that for each galaxy, the [3.6]- 
[4.5] colour is nearly constant with radius and consistent with 
stellar photospheric emission. The 5.8/8.0 ratio (not shown in 
Fig. [3]) is also approximately constant (~0.5), si milarly to other 
galaxies such as NGC 300 bv llHelou et all (2004), even though it 
might not be a general feature. In both the 5.8 and 8.0 /mi bands 
there is an increase of the diffuse emission relative to the 3.6 
and 4.5 /mi bands, as in our Galactic plane. The 8.0/24 colour 
(not shown) also does not show any evident radial gradient or 
prominent structure and is practically constant across the disk. 
Similarly to the 5.8/8.0 ratio, this result points out the overall 
constancy in the ratio of PAHs and warm dust emission in this 
galaxy. 

On the other hand, the 4.5/8.0 images (right panel of Fig. [3]) 
and the 3.6/5.8 ratio (not shown) both retain the structural pat- 
tern of M 33, tracing the inner features, the two main spiral arms 
and the filamentary structure of the disk. In the first case, this is 
probably due to the predominant PAH features in the spiral arms 
and filaments; in the second, the ISM as reflected in the 5.8 /mi 
band is emerging relative to the stars at 3.6 /mi. The 8.0 /mi 
emission along the spiral arms and along filaments surrounding 
the dimmer star forming regions at large galactocentric distances 
become prominent. 

We also constructed an IRAC three colour image of M33 
by subtracting the stellar component in the shorter- wavelength 
images from the 5.8 and 8.0 /mi ones. To remove the stellar con- 
tri bution from the 5.8 and 8.0 ^m imag e s, we used the recip e 
bylHelou et all d2004b : lPahre et al.1 (l2004l) : ICalzetti et al.l (l2005h : 
iPerez-Gonzalez et al.l (120061) . The 3.6 and 4.5 /mi images were 
combined assuming colours appropriate for an M0 III star ([3.6]- 
[4.5]=-0.15 in Vega mag). The final combination is shown in 
Figure0]and displays the stars in blue, the star- subtracted 5.8 //m 
in green and the star- subtracted 8.0 /mi in red. The dominance of 
the 8 /mi emission in the spiral arms is clearly revealed by Fig.ffl 
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Fig. 3. Colour images representing the ratio of the 3.6 over 4.5 jum images (left) and 4.5 over 8.0 jim (right). These images have a 
common resolution roughly equivalent to that in the 8.0 yum image. White corresponds to low colour ratios, and black to high ones. 




The underlying stellar disk is also evident; the excess of stellar 
light in the centre has been pointed out by previous near-infrared 
ground observations dRegan & V ogel 1994), and is related to the 
possible presence of a bulge or bar component. 



3.4. Distribution of cool dust 

Over the last fifteen years, many studies have been devoted to the 
FIR emission of M 33 using IRAS and ISO. iRice et al.l dl990b 
observed M33 with IRAS at 12, 25, 60, 100 jim and found a 
striking correlation between the spatial IR structure and the Hn 
regions. They estimated that 75% of the integrated IR emission 
arises from stars younger than a few times 10 8 yr. The close 
correspondence between the Ha image and the FIR (60 an d 
100 jjm) morphology was also noted by Devereu x et al .(1997). 
iHippelein etafl (l2003h studied M 33 with ISOPHOT at 60, 100 
and 170 /mi, and concluded that the FIR emission is composed of 
a superposition of two components: i) a warm one, found in spi- 
ral arms and SF regions, heated by UV radiation from OB stars, 
contributed by dust at about 45 K; and ii) a cold one around 16 K 
which is smoothly distributed over the disk and heated by diffuse 
interstellar radiation. 

More recently. Ikinz et al.l ([2004) used MIPS observations of 
M 33 to investigate the origin of the FIR emission, and found 
that the emission in the 24 and 70 jim MIPS images follows 
closely the structure of the ionised gas, suggesting heating by 
hot, ionising stars. The overlays of Fig.[6](see Sect. 13.51) show the 
link between 24 and Ha, i.e. between the hot d ust emission and 
the f ormation sites of massive stars (see also Tabata baei et al.l 
2007a). This link is also evident by comparing the Ha and 70 /mi 
emission maps. Moreover, comparing thermal and non-thermal 
radio images of M 33 with MIPS suggests that the 24 and 70 jjm 
emission is more closely associated wi th the thermal radio com - 
ponent than with the non-thermal one (iTabatabaei et al.l l2007a). 



Fig. 4. M 33 three colour image, composite of the central region 
of the complete IRAC data. Blue is a weighted mean of 3.6 and 
4.5 yum, representing the evolved stellar population. Green de- 
picts the 5.8 /mi light after removal of stellar emission, and red 
traces the 8.0 jim channel also corrected for stellar emission [See 
the electronic edition of the Journal for a colour version of this 
figure.]. 



Both the 24 and the 70 jim emission originates in fact from warm 
dust associated with massive star formation. We discuss in more 
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Fig. 5. Ha image of M 33 with the over-plot of 8.0 fim contours. 



detail the connection between 24 fim discrete sources and Ha 
knots below. 

A more diffuse component emission at 160 jim is evident 
when comparing it with the 24 jim map. Cold dust heated by a 
diffuse interstellar radiation field clearly contributes to the longer 
wavelength MIPS band. 



3.5. Comparison with Ha 

In Fig. [5] and Fig.O we show the 8.0 and 24 jim contours on the 
Ha image. Other than the two brightest Hn regions (NGC 604 
and IC 133), the highest levels of 24 fim emission are concen- 
trated in the inner part of the galaxy, and the 24 fim emission 
mainly follows the spiral arms. Notably, not all the Ha bright 
spots have bright 24 jim counterparts. Moreover, the centre of 
the 24 jim emitting sites can be shifted with respect to the centre 
of the Ha emission sites. This might indicate that some of the 
optical emission from the Hn regions is obscured by dust or may 
be due to different intrinsic structures related to the evolution of 
Hn regions. Of particular interest are the 24 jim emission spots 
not associated with Ha emission, which could trace deeply em- 
bedded Hn regions, or a different population of objects. In Fig.[5j 
we can see the same general trends with a concentration of the 
8 yum emission in the centre and along the spiral arms of the 
galaxy. At larger radii, the 8 fim emission becomes more diffuse 
respect to the 24 jim component which matches better the bright 
Ha knots. At 8 yum, the southern arm shows a much more con- 
centrated emission than the northern arm. A similar trend has 
also been found for the molecular gas: giant molecular clouds 
and diffuse molecular gas emission follow the pattern of the arm 
in the south as traced by Ha much more c losely than in the north 
(lEngargiola et al.ll2003l:lHever et al.ll2004l) . 
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Fig. 6. Ha image of M 33 with the over-plot of 24 /dm contours. 

4. The IR emission of nebulae in M33 

We used six catalogues of three classes of known-type sources in 
M 33: three catalogues of Hn regions selected at different wave- 
lengths, two catalogues of SNRs, and one catalogue of PNe: 

1. IR Hn regions: 28 sourc es iGiveon et al . (2002); 

IWillner & Nelson-Patell (l2002h 

2. Radio Hn regions: 40 sources iGordon et all (fl999|) 

3. Optical Hn regions: 78 sources iMagrini et al.l d20 07b) 

4. Optically selected SNRs: 98 sources IGordon et al J dl998b 

5. Radio selecte d SNR s: 30 sources iDuric et aiT dl995): 
IGordon etafl (1 19991) 

6. PNe: 152 sources ICiardullo etafl d2004 

The properties of these will be compared to our catalogue of 
24 yum-selected sources in Sect. [6] 

We have selected our sample of Hn regions from three dif- 
ferent catalogues, which include Hn regions of different mor- 
phologies, locations within the galaxy, and brightness. Our 
IR sample was obtained by m ergi ng the source s obser ved by 
IWillner & Nelson-Patell (120021) and IGiveon et al.l (120021) . It in- 
cludes in total 28 sources located in three chains, to the east, 
west, and south of the galaxy nucleus. Our radio sample com- 
prises 40 Hn re gions with 20-6 cm spectral index (y a ) a > -0.2 as 
catalogued by lGordon et al.l (fl999), and extends out to 20' galac- 
tocentric radii. A list of 78 optically selected Hn regions, from 
narrow band images and with k nown optical spectra, was kindly 
provided to us by L. Magrini (Masrini et al. 2007b). Such re- 
gions are in general much fainter than the IR or radio selected re- 
gions and are located preferentially at large galactocentric radii. 

To probe dust emission during the final stages of stellar evo- 
lution, we selected two samples of SNRs, one de fined optically 
and one defined in the radio. The optical sample (IGordon et al.l 
1 19981) contains 98 SNR candidates selected primarily on the ba- 
sis of optical emission-line ratios. The radio samp le of SNRs 
contai ns 30 (radio-bright) sourc es appearing both in lDuric et al.l 
(119951) and IGordon et al.1 (1 19991) . Their selection criteria are ef- 
fective at finding middle-aged remnants. In fact, since their 
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SNRs have diameters in the range 20 - 48 pc, most of the SNRs 
should be in the Sedov-Taylor expansion stage. The contamina- 
tion in their sample is estimated to be less than 1 , but there are 
20 SNRs embedded within Hn regions. Still, several remnants 
are associated with nearby Hn regions and some are embedded 
in them. 

Fin ally, we adopt the catalogued PNe of ICiardullo et al.l 
(2004). As a result of a photometric and spectroscopic survey, 
they identified 152 P Ne candidates in M 33. By comparison with 
Mag rini et al.l (1200 ll) . the contamination by Hn regions or SNRs 
is estimated at -22%. 

4.1. Photometric methodology 

We have six source lists of known type, for a total of more than 
400 objects distributed over the entire disk of M 33. Our aim is to 
use the IRAC/MIPS-24 colours of these as empirical diagnostics 
to infer the dust properties in various kinds of nebulae. This is 
done by comparing their IR colours with the theoretical colour 
diagrams (see Sect. 0) . 

To these source lists we then add the 24-yum selected sam- 
ple and we shall infer the nature of these sources by compar- 
ing their IR colours with those of known-type of sources (see 
Section [6]). We are also interested in the Ha emission of our ob- 
jects, in particular for the 24-yum selected sample and therefore 
we need to perform photometry on the Ha image. Therefore, 
we must perform photometry in a total of six bands (4 IRAC, 
MIPS-24, and Ha). We merge the IR photometry, and consider 
only those sources which have measurable emission in all the IR 
bands examined, in order to place them on our diagnostic colour 
diagrams. 

The proximity of M 33 combined with the superb resolu- 
tion of Spitzer makes photometry of individual sources in the 
galaxy challenging. Isolating sources is difficult, because of the 
extremely crowded fields, especially the inner disk and circum- 
nuclear regions. Moreover, we expect that many of the objects 
(e.g., Hn regions and SNRs) will be extended, so that we can- 
not accurately perform point- source fitting and subtraction with 
point-response functions. In addition, many classes of sources 
are not prodigious emitters in the mid-infrared (MIR), which 
further complicates the automatization of the photometric proce- 
dures. Hence, we investigated many photometric schemes before 
finally arriving at a reliable solution . 

Following Calzett i et al . (2005), we first attempted aperture 
photometry of the sources in our lists. The initial centres were 
taken from the source lists described in Sect. 01 and the centre 
of the virtual aperture was then allowed to vary up to a certain 
distance from the initial position. While this is a viable solution 
in other situations, it turned out not to be feasible in M 33. First, 
the definition of the fiducial aperture should depend on source 
extent, but this differs among the source classes, so it was not 
clear how to define the limiting aperture. This was an especially 
vexing problem for the wavelengths in which a given source does 
not have measurable emission, since the automatic centring rou- 
tines did not converge. Second, the determination of the back- 
ground was extremely problematic because of the variable de- 
gree of crowding and the variable degree of diffuse emission. 
It was also unclear how to define either "regional" or "global" 
backgrounds in the case of M 33. Thus, we were forced to find 
another solution. 

We then atte mpted elliptical cro ss-section Gaussian fitting 
of sources (e.g. JCalzetti et al.lll995l) . The fitting box was fixed 
to a series of dimensions (from 25 to 10 pixels), and the cen- 
tre, the orientation, the width, and a constant background level 



were left free to vary. As before, the initial centres were taken 
from the source lists. The larger fitting boxes allowed the cen- 
tre to wander too much from the initial position. Hence, the 
source identified by the fitting routine frequently did not cor- 
respond to the source in the list, although reasonably accurate 
backgrounds were obtained when checked against the source 
surroundings. The smaller fitting boxes better constrained the 
centre, but in crowded fields gave a background which was much 
too high. This resulted in gross underestimations of the total flux, 
especially when the source was rather faint and comparable to 
the estimated background. An additional problem encountered 
with Gaussian fitting was its unsuitability for multiple sources. 
Many of the sources in our lists are in truth "clusters" of smaller 
sources when examined with the Spitzer spatial resolution and 
at wavelengths different from those in which the sources were 
defined. These problems taken together led us to distrust this 
method for a large fraction of sources in our lists. 

To overcome these shortcomings, we finally devised a "hy- 
brid" solution, which optimally combines the photometry from 
Gaussian fitting and that from isopho tal photometry. This last 
was accomplished with SExtractor (Ber tin & Arnoutsl Il996l) 
which is able to deal with very large images automatically, and 
handles reasonably accurately a wide variety of object shapes 
and sizes. The idea was to calculate the total flux integrated 
over an irregularly shaped aperture containing flux levels above 
some isophotal limit. The hybrid method consists of three steps: 
(1) perform Gaussian fitting of all sources as described in the 
previous paragraph; (2) perform SExtractor photometry as de- 
scribed below; (3) merge the two sets of photometry according 
to an algorithm which decides which photometric measurement 
to prefer, in order to have one measurement per source. 

For (1), we adopted a fitting box of 10 pixels; this was 
the best compromise between centre wandering and background 
level, as the photometry was typically centred on a position rel- 
atively close to the initial one, although with a potentially spuri- 
ously low flux because of the background determination. 

In practice, for (2), we used the ASSOCIATIVE feature 
of SExtractor, and input our source lists with initial positions. 
The initial SExtractor photometry consisted in those sources ex- 
tracted at a 1 Oct level above the background, at a distance within 
15 pixels (= Ar) from the initial position. These two values re- 
sulted to be the best compromise between source association 
and source identification. Lower values of signal-to-noise ratio 
(< 10cr) resulted in more sources being identified, but were in 
the mean farther away. The background for the SExtractor pho- 
tometry was defined to be the local median of the image in a 
64-pixel box centred on the source. 

Step (3), the merging of the Gaussian fitting and SExtractor 
photometry, was accomplished by examining the radial dis- 
tance of each measurement from the initial position. When the 
SExtractor centre was within Ar of 5 pixels (in radius), and as 
long as it was within a factor of two of the Gaussian fitting cen- 
tre, it was preferred over the Gaussian method. We experimented 
with a few values of Ar and a few values of the relative ratio of 
the radial distances, and concluded that this gave the best results 
overall. Indeed, it appears that this hybrid approach overcomes 
the problems of underestimating source fluxes because of back- 
ground estimates, and successfully calculates the total flux even 
when sources are in fact multiple. Indeed, the SExtractor pho- 
tometry and the Gaussian fitting results are fairly well correlated 
except for low fluxes with high background where Gaussian fit- 
ting tends to underestimate the true flux, and high fluxes for ex- 
tended or clustered sources. 
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Table 1. Detection rate in all four IRAC and IRAC+MIPS-24 
bands. 
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Fig. 7. Detection rate in the different wavebands. Open symbols 
represent the Hn regions selected in IR (squares), radio (trian- 
gles) and optical (circles). The (blue) crosses represent the opti- 
cally selected SNRs and the (blue) plus signs the radio selected 
SNRs . The filled (red) diamonds represent the PNe. The (green) 
tripods represent our 24 jim sources. 



4.2. Detection rates 

The detection statistics are illustrated in Fig. [7] and Table [T] The 
matching radius Ar max for the coincidence of sources, was deter- 
mined by examining the scatter on the diagnostic diagrams, see 
Sect. [5] In the end we adopted a value Ar max = 3 pixels (about 14 
pc), which is sufficiently large to accommodate possible astro- 
metric uncertainties. 

The positions of catalogued sources with photometry in all 
four IRAC bands are shown in Fig.[8j on the IRAC 4.5 jim image 
of M 33. The upper panel refers to Hn regions and we note that 
the optically selected sample is mainly in the outer regions of 
the galaxy, while the IR and radio selected Hn regions are in 
the inner parts of M 33. The lower panel shows the locations of 
PNe and SNRs. The single detection of the radio- selected SNR 
coincides with an optically- selected one. 



4.2.1 . Hn regions 

The detection rate of Hn regions depends on the waveband and 
on the catalogue. As expected, the brighter Hn regions (radio and 
IR sample) are detected at a higher rate than the faint ones (op- 
tically selected catalogue). The detection rate for bright regions 
increases even more if we loosen the strict requirement about 
the centre displacement. For large Hn regions in fact the centre 
of the IR emission can be displaced by more than 3 pixels from 



the nominal position, and hence not accepted as coincident. The 
lower resolution in the 24 jim band makes objects rounder, and 
the Gaussian fitting is often more accurate than in the 5.8 and 
8.0 /mi band. This partly accounts for the higher detection rate 
in the 24 jim MIPS wavelength. But, there is a general tendency 
for the Hn regions to be less visible in the shortest wavelengths 
(3.6 and 4.5 jim) than in the 8.0 and 24.0 jim channels. The de- 
tection rate at 5.8 //m always shows a drop, with respect to the 
adjacent channels (see Fig. [7]). 

4.2.2. Supernovae remnants 

We found a low rate of IR counterparts for the SNRs identified 
in optical and in radio. The percentage of detections is about 
15% for both samples. Four SNRs were detected in at least one 
band: two were detected at all wavelengths, and other two only at 
4.5 jim. The four detected SNRs have radii of ~ 8 - 16 pc, while 
the two undetected ones have radii greater than 20 pc and may 
be of a lower surface brightness. The low number of SNR de- 
tections is consistent with previous IRAS st udies trying to find 
Galac tic SNR counterp arts in the IR: 17% (lArendtjll989h and 
18% (ISaken etalJI 1992b . Also, using IRAC. lReach etal.1 (120061) 
found only 18 IR counterparts out of 95 SNRs in the inner 
Galaxy. Very few Galactic S NRs were detected by MSX, the 
Midcourse Space E xperiment (iCo hen et al. 2006). In the Large 
Magellanic Cloud, feorkows ki et al.l (120061) f ound no emission 
from Type la SNRs and lWilliams et aLH l2006) did not detect the 
two largest SNRs of their sample in any of the IRAC or MIPS 
wavebands. 

In our samples, the detection rate is also about 30% less 
in the 8.0 and 24 yum bands with respect to the shorter IRAC 
wavelengths. This might be due to the destructive effect of the 
blast wave on the smaller grains; PAHs can already be razed 
in slow shocks and they a re never detected in dense, shocked 
clumps (iReach et al.ll2002l) . The emission seen at 3.6 and 4.5 jim 
is likely due to ionic and molecular line s, such as Bra 4.1 fim 
(Re ach et al.l 120061: IWilliams et al.l 120061) . The environment of 
the SNRs should have a strong influence on their IR detection 
rate: remnants interacting with a denser external medium should 
be more easily detectable. If PAHs in the surrounding IS M are 
destroyed by nearby supernova events (Gor jian et al .1120041) . they 
will not emit at 8.0 yum (nor contribute to 3.6 and 5.8 yum), which 
could explain the quite low rate of detection of a remnant in all 
4 bands. 



4.2.3. Planetary nebulae 

The detection rate of PNe is about 25% for the three shorter 
IRAC wavelengths and then rises steeply to reach more than 
50% in the 24 /im band . This is consistent with the result by 
ICohen & Barlowl (120051) who detect PAH emission i n less than 
half o f a sample of 43 Galactic PNe. In addition, iHora et all 
( 2004) claim that the IRAC colours of PNe are red; despite typ- 
ically negligible continua, they are bright in the 8.0 jim band 
probably because of a sizable contribution from the two H2 lines 
and an [Ar III] line at 8.99 jim. 



4.3. Spectral energy distributions 

The average IR Spectral Energy Distribution (SED) for our dif- 
ferent lists of objects is shown in Fig. [9j Hn regions, either IR 
or radio selected, display the highest fluxes, about two order of 
magnitudes above the faintest in the list, that is the PNe. The 
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Hn optical sample instead, largely consisting of quite faint re- 
gions, sits at lower values, in between the two samples of SNRs. 
The steep rise between 8.0 and 24 /mi of all the Hn regions sam- 
ples, implies that Hn regions are invariably associated with warm 
dust. SNRs and PNe show milder slopes, their SEDs being al- 
most flat, or even declining as for optically selected SNRs. This 
last feature reveals changes of the ISM characteristics related to 
SN events. 

In Hn regions, the drop between 3.6 and 4.5 fim can be 
accounted for by the photospheric contribution of young stars 
in the two shortest wavelength IRAC channels or perhaps line 
emission. Candidate lines in the 3.6 /urn band includ e [Fe m] 
3.229 yum, and [Co m] 3.492 /zm (iGerardv et al.ll2007l) . but line 
emission can be important in other bands as well as long as there 
is little continuum e m ission. This drop has a lso been detected b y 
Mer cer et al.l (2004): Goriian et al. (2004k Ijones etal.1 (120051) : 
Williams ~al.l (l2006l) : lGerardv et al.1 (120071) . The absence of the 
4.5 jim drop in the means of radio selected Hn's is due primarily 
to a single object with a high 4.5 /mi flux. 

5. Diagnostic diagrams 

5.1. Theoretical diagnostic diagrams for IRAC and 
IRAC/MIPS flux ratios 

To better understand the IRAC/MIPS colour diagrams of dis- 
crete sources, we assembled a series of models with different 
physical parameters, as described more in detail below. The IR 
colours are derived by convolving the model spectra with the 
IRAC and MIPS-24 filter response curves, as provided by the 
Spitzer Science Center web page. The colours are here computed 
as ratio of flux densities (in Jy, i.e. per unit frequency). 

5.1 .1 . Hn regions with Cloudy 

Hn regions are possibly the most prominent class of sources in 
our images. We modelled them with version 5.07 of Cloudy , 
last described bv iFerland etaD (1 1998b (see also lFerlandll2003l) 
The c entral star was taken to be a log g=4 Kurucz model (iKuruczl 
1979) at T e ff = 37 000 K, as modified by Rubin to better account 
for the stellar flux in the model atmospheres between 41 eV and 
54 eV. In all cases, we assumed a stellar ionising luminosity of 
10 49 sec" 1 , spherical geometry, constant gas density, a covering 
factor of 0.005, and a metal abundance of 50% solar. 

Grains were "turned on" in the Cloudy models so that heat- 
ing and cooling of the grains and the gas are calculated in a 
self-consistent way. A standard PAH distribution is included in 
the grain population, with tr eatment of photoe lectric heating and 
collisional processes (see IFerland et al.lll998l) . We did not ex- 
plicitly introduce a photo-dissociation region (PDR) in the cal- 
culations. 

We accounted for the range in typical Hn region parameters 
by varying the outer radius and the total hydrogen density so as 
to roughly follow th e well-known size-den s ity n e ~ D~ l rela - 
tion for Hn regions (Garav & Lizano 1999; Kim & Koo 2001). 
Four models were calculated with the following pairs in size- 
density space with density given first in cm -3 followed by outer 
radius in pc: (1000,2); (100,10); (20,30); (10,60). 

The resulting spectra are shown in Fig.[l0l where the (blue) 
short-dashed line shows (100,10); the (cyan) dotted line (20,30); 
the (green) long-dashed line (10,60); and the (red) solid line 
(1000,2). It can be seen from the figure that varying the den- 
sity by a factor of 100 drastically changes the PAH emission and 
its dominance relative to the continuum. Since the PAHs in most 



Hn regions originate in t he surrounding PDR rather than in the 
ionised region itself (e.g I Whitney et al .1120041; IChurchwell et al] 
120041; LDopita et al.ll2006l) . the simulations may not be truly rep- 
resentative of Hn regions unless a PDR is absent. However, 
they clearly span a wide range in IR peak wavelength and 
PAH strength, so should be adequate for helping to interpret the 
IRAC/MIPS24 diagnostic diagrams. 

5.1 .2. Embedded Hn regions with DUSTY 

The Hn regions we have modelled with Cloudy are basically 
optically thin, because we have not performed multiple iter- 
ations on the solution. Moreover, even with iteration, optical 
depth effects in the dust itself are not taken into account in the 
Cloudy formalism. Hence, to better model compact and ultra- 
compact Hn regions which are expected to be embedded in an 
opaque dust cocoon, we have run a series of models with DUSTY 
ivezic &Elitzur 1997). DUSTY solves semi-analytically the ra- 
diative transfer physics in a dust envelope by exploiting the self- 
similar nature of the problem. The dust is located in a shell ex- 
ternal to the Hn region itself. The program allows for a variety 
of spherical shell configurations, and provides a realistic spread 
of dust temperatures as the radiation field changes with distance 
from the source. At the same time, this code suffers from some 
limiting assumptions, such as spherical symmetry, no treatment 
of stochastic processes for PAH emission, and the decoupling of 
the gas from dust. Nevertheless, for our purposes, it is a power- 
ful tool with which to investigate the dust continuum emission 
from young dust-enshrouded stars and star clusters. 

With DUSTY, it is necessary to specify the temperature at 
the inner radius of the dust shell, Ti n which, for a given incident 
radiation field, determines the inner radius of the region, R{ n . The 
radial distribution of the dust is described through the ratio of the 
outer radius to the inner one, Y out (= R ut/Rm), and the power- 
law exponent p describing the radial fall-off RP . In all cases, 
we a dopted a standard "MRN" size dist ribution dMathis et al.l 
11977b . and a standard dust composition (Draine & Led 119841) . 
Finally, we must specify the radial optical depth Ay which, for 
a given geometry, fixes the dust mass. The DUSTY models were 
calculated for 10 values of Ay, ranging logarithmically from 0.1 
to 100 (specifically 0.1, 0.215, 0.464 for each decade in Ay). We 
adopted single-zone models, with only one stratum in the shell. 

Two temperatures were input for Ti n , 300 K and 600 K, 
which correspo nd roughly to th e range observed in extragalac- 
tic Hn regions (iHunt et al] |2005). Because of the non-negligible 
thickness of the dust shell, such relatively high temperatures 
at its inner rim are also required to explain the average tem- 
peratures observed in G alactic Hn regions (iConti & Crowtherl 
2004l: lGiveon et al.ll2007l) . We fixed the power-law exponent p to 
-0.2 as a compromise between uniform regions and those with 
a steeper fall-off. This exponent strongly affects the amount of 
cool dust in the shell, but does not significantly influence the 
MIR portion of the spectrum that we are examining with the 
IRAC/MIPS-24 diagnostic diagrams. The thickness ratio Y ou t 
varied from 10 to 10000; Y ou t also controls the amount of cool 
dust in the shell as thicker shells have proportionately more cool 
dust. 

To model the exciting star cluster input to D USTY, we in- 
corp orated the Starburst 99 simulations (SB99, iLeitherer et al.l 
1999). Because we are interested in young, dusty star clusters, a 
single age of 3 Myr was used. SB99 takes into account the neb- 
ular continuum emission which can dominate the near-infrared 
(2 fim - 4 //m, NIR) spectral range, and thus could in principle 
significantly influence our IRAC colour diagrams (see below). 
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In practice, the slope of the input spectrum of a 3 Myr starburst 
is quite similar to the star used for the Cloudy simulations, at 
least down to 0.1 yum, so we expect few systematic differences 
in the dusty star clusters modelled with DUSTY relative to the 
Cloudy models. The metallicity of the input spectrum was taken 
to be 20% solar; this was to better account for possible lower- 
metallicity embedded sources in the outer disk. However, as il- 
lustrated below, this choice of metallicity does not significantly 
affect the resulting IRAC/MIPS-24 colours. 

We show examples of the resulting DUSTY spectra in Fig. 
ITTI Here we plot only 5 values of Ay, namely 0.1, 2.15, 10.0, 
21.5, and 100 (shown as long-dashed, short-dashed, dotted, dot- 
dashed, and solid lines, respectively). Also shown (as a heavy 
solid line which ends at X ~ lOOyurn) is the input stellar spec- 
trum of a 3 Myr starburst from SB 99. The figure clearly shows 
the importance of nebular continuum emission and the onset of 
warm dust in the spectral region sampled by IRAC. 

5.1.3. Blackbodies 

To illustrate the similarity between the emission from stellar 
populations and blackbodies in the IRAC bands, we also calcu- 
lated blackbody spectra at temperatures ranging from 100 000 K 
to 300 K. In the IRAC/MIPS spectral region, there is virtually no 
difference among the blackbodies from 100 000 K to 3500 K. As 
expected, their colours are also very similar to the SB 99 simula- 
tions of the stellar populations in a 3 Myr cluster. 

5.1 .4. Interstellar medium with PAH features 

An example of the spectrum of the diffuse ISM at high Galactic 
latitudes was also includ ed in the diagnostic s. The spectrum was 
taken from Figure 8 of|LL&Draine (200lb. and the digital ver- 
sion was kindly provided by B. Draine. In the IRAC/MIPS-24 
spectral range, the diffuse high-Galactic latitude ISM shows ex- 
tremely strong PAH emission, superimposed on a rising contin- 
uum. 

5.1.5. Reflection nebula NGC7023 

NGC 7023 is a well-studied PDR that, because of its proxim- 
ity, is an excellent source to study the physical processes occur- 
ring in such regions. It is a reflection nebula illuminated by a 
B3Ve star, and is located in a cavity of the molecular cloud. This 
cavity ho sts a dense PDR w ith its peak located away from the 
star itself dFuente et al.l2Q0Qh . We obtained the short- wavelength 
portion of the ISO s pectrum of this ob ject from the uniform 
database published bv lSloan et all (l2003l) . Like the diffuse ISM, 
NGC 7023 also shows relatively strong PAH features superim- 
posed on a rising continuum. 

5.1.6. Results 

The IRAC and IRAC/MIPS-24 diagnostic diagrams relative to 
the models described above are shown in the bottom panels of 
Fig-El Fig- [13 and Fig. [I4j Blackbodies are depicted by filled 
black circles. NGC 7023 is represented by a (yellow) triangle, 
the diffuse high-latitude ISM by an inverted (cyan) triangle, the 
SB 99 simulations by (grey) diamonds. The Cloudy Hn region is 
represented by filled green squares. Finally, the DUSTY models 
are represented by filled (Tin = 300 K) and empty (Ti n = 600 
K) circles, for a variety of shell thickness Y out = 10 (solid), 100 
(dashed), 1000 (dash-dot) and 10000 (dotted line). 



There are several interesting features that can be noticed 
in the plots. First, a PAH-dominated spectrum such as the dif- 
fuse ISM, the small dense Cloudy Hn region, and the reflec- 
tion nebula have virtually identical IRAC colours to some of the 
Ti n =300 K DUSTY models. The implication here is that strong 
PAH features can be easily confused with a rising continuum 
in the MIR. The only way to resolve the ambiguity is through 
detailed spectral information or by using the 24 \im filter in con- 
junction with the IRAC 8 [im one. Indeed, our plots show that 
only 2 IRAC filters, together with MIPS-24, are sufficient to 
separate strong continuum emission in the MIR from strong 
PAH features. Second, Figure [121 shows that strongly embed- 
ded sources (e.g., Ay ~ 100) occupy a well-defined region in the 
IRAC colour diagram. This may potentially be useful for sin- 
gling out ultra-compact dusty Hn regions from our 24 yum source 
catalogue. Third, stars (roughly equivalent to 3500 K blackbod- 
ies) occupy a well-defined region in both diagrams, and should 
therefore be easily separable from other classes of objects. 

On the negative side, the DUSTY model predictions in the 
IRAC/MIPS-24 diagram are not easily distinguishable from 
Cloudy predictions. Moreover, even the highly embedded clus- 
ters do not occupy a distinct region in the colour diagrams. This 
makes the use of IRAC bands necessary, in order to better sepa- 
rate high extinction regions from optically thin ones. 



5.2. Observational diagnostic diagrams 

A comparison of the IRAC colours of the models + templates vs. 
observations of the different kinds of nebulae of Sect. [His shown 
in Figs. [12l and \13\ Diagnostic diag rams using IRAC colours 
were presented by lReach et al.l (120061) (focusi ng on SNRs emis- 
sion mechanisms) and by ICohen et al.l (l2007|) . We present two 
versions of the IRAC diagram, in order to better assess the di- 
agnostic value of the different ratios. We also propose in Fig.[T4l 
a new diagnostic diagram involving the 24 fim MIPS channel 
which allows us to disentangle some degeneracies present when 
using only IRAC colours. In the next Section, will use both 
the theoretical colour diagnostics and the observed colours of 
known-type objects to investigate the nature of sources selected 
at 24 /im. 



5.2.1 . Hn regions 

In the first IRAC colour diagram (see the middle panels of 
Figs. [TJland [13]), we find a concentra tion of Hn region s in the 
zone dominated by PAH emission (see lReach et al.ll20 06). In the 
same zone, lie the reflection nebulae NGC 7023 and the diffuse 
Galactic ISM as well. The dispersion in the colours of Hn re- 
gions makes it difficult to separate them from SNRs and PNe. 
Statistically, the Hn regions are more concentrated near the re- 
gions defined by low values of the two shortest IRAC channels 
and/or equivalently high values of the two longer IRAC wave- 
lengths. The most intense PAH bands (6.2, 7.7, 8.6 jim) are lo- 
cated in the 5.8 and 8.0 jim channels and this may explain the 
colours of Hn regions. 

We compared the IRAC colours of ou r catalogued Hn re - 
gions with the diffuse Hn regions studied b vlCohen et al] (120061) . 
Upon conversion to their magnitude scale dPahre et al.ll2004|) . we 
find the IRAC-colour distributions to be mutually consistent. 
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5.2.2. Supernovae remnants 

As a word of caution, since (especially young) SNRs are often 
associated with Hn regions, they are often hard to discern and 
to measure photometrically given the crowding and confusion 
in the field. In the colour diagrams, the few detected SNRs lie 
between the Hn regions and the PNe, generally closer to the Hn 
regions. 

Since the colours of pure synchrotron emission are never ob- 
served, there are several alternative emi ssion mechanisms f or the 
IR radiation that have been proposed by lReach et al.l (2006) such 
as shocked molecular lines, shocked ionic lines, and PDR PAH 
emission. About half of their SNRs have colours consistent with 
molecular shocks, 17% are consistent with ionic shocks, and 
22% with PAH emission from unshocked ISM. The others are 
thought to represent a mixture of shock types. All our 4 SNRs 
appear to belong to the class of shocked molecular emission (see 
the -grey- filled rectangle in the middle panel of Fig. [12]), that is 
shocked molecular gas cooling mainly via MIR emission lines. 
Also, optically, SNR s appear as shock-heated emission nebulae 
dGordon et al.l 1999k . con firming the p r evalen ce of the two shock 
mechanisms described by lReach et al.l (120061) . The colours of our 
SNRs are also close to the colours of some Hn regions. Note that 
the colours of the SNRs could depend of the epoch at which the 
SNR is observed (for instance because of the onset or decrease 
of shocks as they expand into the ISM). 



5.2.3. Planetary nebulae 

With few exceptions, catalogued PNe are confined to a region 
characterised by low values of the 5.8 and 8.0 yum emission, the 
bandpasses where PAHs are believed to contribute substantially 
to the IR emission. Indeed, PNe spectra usuall y show an absence 
or very weak presence of P AH, as noted by Bern ard Salas et al.l 
(200 2b with ISO and bv [Bernard-Salas et al.l (120041) and with 
Ber nard-Salas et al.l d2006h Spitzer. In the 4.5 ^m band, the [Mg 
iv] 4.485 yum and [A vi] 4.530 yum lines could give an ad ditional 
contribution for high-excitation PNe (iCohen et al. 1 12005). Thus, 
in absence of detailed spectral information Spitzer colours can 
prov ide some diagnos tics to help to identify PNe candidates (see 
also iHora et al.ll2004 . 



6. The 24 yum sources in M 33 

Spitzer images provide a unique opportunity to compile vast cat- 
alogues of MIR discrete sources down to stellar luminosities for 
galaxies in the Local Group. In the present section we describe 
and discuss some features of the catalogue of 24 yum sources 
(24Ss) which we have generated for M33. In the first subsec- 
tion, we describe the source selection criterion and flux measure- 
ments; in the second subsection, we investigate their nature in 
light of both theoretical and observational diagnostic diagrams; 
in the last subsection we present their luminosity function. 



515 objects.The brightest sources at 24 jim are located in the 
central region and along the spiral arms. Fainter sources are more 
smoothly distributed but still follow the flocculent pattern of the 
galaxy. 



6.2. Nature of the 24 jim sources 

Perhaps not surprisingly, sources in our 24 //m catalogue turn 
out to have a higher detection rate than any of the other kinds of 
objects considered in the previous sections. This is particularly 
marked at 5.8 yum where the detection rate for the 24 /mi sample 
is much higher than for Hn regions (see Fig. [7]) although, as seen 
in Fig. [9j the averaged SED of the 24Ss is similar to the one of 
Hn regions. 

The position of the 24Ss in the colour diagrams are illus- 
trated in the top panels of Figs. [12 [13] & H In Fig. [121 the 
region with the highest concentration matches well the one of 
Hn regions. However the distribution of 24 /mi sources extends 
into the regions populated by SNRs and PNe. The comparison 
with the theoretical models indicates a paucity of sources with 
high extinction. However, as Fig. [12] shows, neither a pure black- 
body, Cloudy, or the two DUSTY models can alone reproduce 
the colour distribution of the 24 yum sources. Fig. [13] shows a 
clear departure of the 24Ss sample from the expected narrow 
4.5/5.8 colour dispersion. Such departure is particularly marked 
in the red tail of the distribution where both 24Ss and known- 
type sources tend to have colours similar to the PDR observed in 
the reflection nebulae NGC 7023. 

The hybrid IRAC-24 fim MIPS colour diagnostics of the 
selected sources are illustrated in Fig. [141 No 24S displays a 
24/8.0 yum flux ratio as high as the DUSTY models with shell in- 
ner temperature Ti n = 300 K; the warmer models (Ti n = 600 K) 
seem more appropriate to describe their physical conditions. The 
range spanned by the 24/8.0 ratio in the 24 yum sources is large 
but similar to the one covered by the catalogued Hn regions, in- 
dicating that in both samples the contribution from PAH fea- 
tures may vary widely. So, 24Ss with both a low 4.5/8.0 and a 
low 24/8.0 ratio could well be Hn regions with high PAH emis- 
sion (see, for instance, the location of the radio- selected Hn re- 
gions). The Cloudy Hn region simulations seem to match well 
the colours of the 24Ss, in spite of their large scatter. On the con- 
trary, the observations rule out most stellar colours: both pure 
blackbodies and the SB99 templates. However, as above, the 
general agreement between any single model and the 24Ss dis- 
tribution in the colour diagrams is unsatisfactory; such a trend 
is probably due to a variable PAH component which is not con- 
sidered in most of our templates (but see loci of NGC 7023 and 
the ISM). In the next section, we will select a subsample of 24Ss 
associated with the Ha emission. These sources, likely to be op- 
tically visible Hn regions, achieve a better agreement with low- 
extinction models. 



6.1. 24 fim source catalogue 

To select and measure individual sour ces from the 24 pm im - 
age, we used the SExtractor software (Ber tin & Arnoutsiri996l) . 
Using a convolution with a Gaussian filter we claim a detec- 
tion when the source area is larger than 10 pixels and the peak 
signal-to-noise ratio is at least 10 times the local background. 
Guided by the galaxy pattern traced by the 24 yum diffuse emis- 
sion, we removed spurious sources such as foreground stars or 
artifacts due to map edge effects. The final catalogue comprises 



6.3. Luminosity function of 24 jim sources 

The 24 yum emission of M 33 within an ellipse of major axis 40' 
and axial ratio of 0.62 is 40.9 Jy. At a distance D = 840 kpc, 
such ellipse corresponds to a circle at galactocentric distance of 
4.8 kpc inclined by 52° with respect to our line of sight. The 
24 yum integr ated emiss i on is in good agreement with the value 
reported by iHinz et all (120041) (see their Fig. 2) and implies a 
24 jam luminosity (L = vF v x AnD 2 ) of 4.3 x 10 41 ergs s" 1 . 
The contribution of the 515 discrete 24Ss amounts to 1.4 x 
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10 41 ergs s , hence the diffuse component accounts for two 
thirds of the 24 jim emission. 

The cumulative luminosity function of the 515 discrete 
sources catalogued at 24 jim is shown in Fig. [T6j In this fig- 
ure, the abscissas are reported as flux density at 24 //m in mJy 
and also as total IR luminosity (TIR). To estimate the TIR lumi- 
nosity, we follow the prescription of ICalzetti et al.l (120051) (see 
their Eq. 1). The conversion involves the 24 yum and the 8 fim 
measurements. Since we do not have a detection at 8 jjm for 
all 24Ss, we compute the average 8 yum/24 jim flux ratio for the 
sources detected in both bands: the average is close to unity and 
therefore, in the lower panel of Fig.[l6l which includes all 24Ss, 
we use logL(TIR) = logL(24 yum)+0.908. As a check, in the up- 
per panel of Fig.[T6j we plot the cumulative luminosity function 
only for the sources detected in both bands, using for each source 
the measured 8 yum/24 jim flux ratio: the shape and meaning (see 
below) of the two curves is quite compatible. The main effect, 
when using a constant average flux ratio, is a slight expansion of 
the range in L(TIR) by about 0.4 dex, together with a brightening 
by the same factor. 

From the flattening of the faint end of the distribution, we 
estimate the completeness limit of the catalogue around 1 mJy, 
or an L(TIR) of 5 x 10 37 ergs s" 1 on the aver aged scale. That 
is the bolometric luminosity of a B1.5 V star (ICoxl l2QQQh and 
so our listing should be quite complete even for faint obscured 
Hn regions. Apart from the levelling off at low fluxes, as of- 
ten observed for Hn reg ions, open clusters, and associations 
(McKee & Williams 1997), the luminosity function displays a 
double slope behaviour, markedly steeper at the high luminosity 
tail. The change of slope, in the simplest scenarios, represents 
the change of regime between poor and rich clusters, where rich 
means numerous eno ugh to represent fairly the high-mass IMF 
(lOey & Cla rke 1998). In this framework, the transition point be- 
tween the two regimes marks the luminosity of the single bright- 
est star; below this value, the observed statistics is modified by 
the sampling variance. We find that the transition point occurs 
around F y (24 fim) = 70 mJy, that is, in the averaged scale, an 
L(TIR) ^ 5 x 10 39 ergs s" 1 which, still b olometrically, corre- 
sponds to an 03 V star (Vacca et al.| [l996). This implies: first, 
that most of the bright 24Ss are in fact luminous young stel- 
lar clusters and, second, that our L(TIR) scale cannot be grossly 
mistaken. From F y (24 jjm) = 1 to F y (24 jjm) = 70 mJy we fit a 
power index -0.55 ± 0.01 (plain -red- line) and, for larger values 
an index -1.09+0.08 (dotted -green- line). It is easily shown that 
an index about unity, in the log-log cumulative distribution, im- 
plies for the rich galactic clusters a distribution N(L)dL oc L~ 2 dL 
or, equivalently, N(n*)dn* oc n~ 2 dn*, n* being the number of stel- 
lar members. Again, this is precisely the slope observed in the 
lumin osity function of Hn regions in late spirals (ICaldwell et al.l 
Il99lh . as well as the slope determined for G alactic clusters and 
associations, e.g. McKee & Williams (1997). 

The brightest source reported in Fig. [161 at F y (24 jjm) 
1.8 Jy, would be roughly equivalent to an open cluster with 
around 3.5 x 10 7 L©; assuming that all members (0.06 to 90 M©) 
are on the main sequence, this makes a cluster of some 10 5 M , 
more than 10 6 members, and more than 1000 ionising (> 15 M©) 
stars. On the other side, the faintest detected discrete 24Ss have 
F y (24 yum) of 0.1 - 0.2 mJy, that is L(TIR) * 10 37 ergs s" 1 . This 
is about the bolometric luminosity of zero age main sequence 
B2 stars (Panagia 19731) : so it is possible that these objects are 
also (the faintest) Hn regions and our luminosity distribution ap- 
pears to precisely match the whole range expected for ionised re- 
gions and their complexes. Anyway at low luminosity, especially 
below 10 4 L , the situation may be more intricate since also 



evolved stars can contribute substantially, more or less depend- 
ing on the history of the region, to the population of IR sources; 
carbon giants, for example, at the distance o f M 33 are expected 
to em it an average F y (24 /mi) ^0.25 mJy dGroenewegen et all 
120071) . The nature of the 24Ss will be further explored, by inves- 
tigating their spatial location within the galaxy and IR colours, 
in a forthcoming paper. 

7. Star Formation 

7.1. Global star formation rates 

In regions of moderate visual extinction, Ha line emission is ex- 
pected to be a reliable indicator of recent star formation. M33 
has globally a low extinction rate and therefore we can assess 
the performance of MIR emission as a star formation tracer, by 
comparing the SFR inferred from MIR emission with that de- 
rived from Ha emission. 

Ha fluxes w ere measured on the map of 

iHoopes & Walterbosl d200Qb cited earlier (Sect. l23l) and 
converted to SFR following Kennicutt (1998) We find a total 
SFR from Ha, within 5 kpc from the centre, of 0.22 M Q yr" 1 
if no extinction corrections are taken into account. On average , 
this is a reasonable assumption since Tabatabaei et al. (2007b) 
found a perfect 3.6 cm-Hd correlation for the thermally emitting 
regions in M 33. In any case, giv en that the diffuse Ha emis sion 
amounts to 40% of the total (IHoopes & Walterbosl kOOOl) . an 
extinction correcti on of individual sourc es b y Ay = 1 on average 
(as c laimed by iDevereux et al.l 1 19971 and llsrael & K ennicutt 
Il980h . would increase the SFR in the same area to 0.3 M© yr" 1 . 

To estimate the SFR from the IR emission, we used two in- 
dependent methods: we first estimate the TIR luminosity within 
a radius of 5 kpc using the data at 8.0 and 24 /mi. With a 
flux of 47.3 Jy, the integrated 8.0 jim luminosity L v (8.0 jjm) 
is 4.0 x 10 28 ergs s" 1 Hz" 1 , corresponding to vL y (8) = 1.5 x 
10 42 ergs s" 1 . This is slightly higher than L y (24 jjm)= 3.5 x 
10 28 ergs s" 1 Hz" 1 (or vL v (24) = 4.4 x 10 41 ergs s" 1 ). Such a 
difference in the 8 and 24yum flux (47.3 vs. 40.9 Jy) is consis- 
tent with the IR SEDs measured in nearby galaxies (iDale et al.l 
2005), in which IR fluxes at 8.0 um often exce e d thos e at 24 fim. 
Using these values, Eq. 1 in ICalzetti et al.l (120051) yields an 
L(TIR) = 1 .Ox 10 9 L Q . The I R SFR is usual l y com puted from the 
FIR (40-120 fim) luminosity. ICalzetti et al.l OQQQh found that the 
ratio Lir( 1-1000 yum)/L IR (40-120 fim) is about 1.75. In our case 
this translates into a FIR luminosity of 5.8 x 10 8 L within a ra- 
dius of 5 kpc in M 33. This is in good agreement with the value 
of 6.5 x 10 8 L found from the total IRAS fluxes reported by 
iRice et all (fl990), especially considering that we are including 
only the central 10 kpc in our estimate. Following the precepts 
of iBuat & Xul (119961) for late spirals, that is SFR [M Q yr" 1 ] = 8 x 
10" 44 L(FIR) with uncertainties of around 50%, we obtain a SFR 
equal to 0.18^ 18 M© yr" 1 . Within the uncertainties this estimate 
is compatible with those derived from Ha emission (with and 
without extinction correc tions). We obtai n the same value from 
the TIR-SFR relation bv iKennicuttl (Il998h . taking into account 
the luminosity conversion (TIR vs. FIR). 

The second method to derive the SFR from the IR emission 
involves the linear relation between radio and IR luminosities, 
and SFR (lYun et al.l200ll) . IWu et all (120051) derived SFR directly 
from the 24 fim luminosity (see their Eq. 1): in our case, this 
translates into a SFR of 0.17 M© yr" 1 . Using the 8 fim lumi- 
nosity (their Eq. 2) gives a higher value of SFR, 0.26 M yr" 1 , 
even after subtracting the stellar contribution to the 8 fim flux, 
-10%. Although within the errors of the 24 fim and FIR results 
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ffrom lRice etaDll990h . the 8 pm value is ^40% higher than ei- 
ther the 24 um or t he Ha one. A similar trend was found by 
ICalzetti et all (120051) in which the 8 pm SFR estimate can be 
significantly larger in weakly ionized regions such as the disk of 
M 33 taken globally. This is because the 8 pm emission is pro- 
duced by more than one mechanism, and thus does not trace SFR 
as accurately as other IR wavelengths. 

In general then both MIR prescriptions for computing the 
SFR reach a satisfactory agreement with the Ha-derived value 
and, beyond confirming the adopted SFR computing scheme, 
implies: a) that actively star forming sites are the main contrib- 
utors to the MIR emission; b) that mean extinction is low also 
in the inner Hn regions; c) that highly obscured sources do not 
contribute significantly to the MIR emission. 

7.2. Star formation rates in 24 pm sources 

We then searched for Ha emission at the location of the 24Ss. 
This was done in order to compare locally the SFRs from Ha 
with those from the MIR, but also to check whether sources 
with and without an Ha counterpart have on average the same 
IR colours. On the narrow-band (continuum subtracted) Ha 
image, we performed the photometry of the 24Ss using the 
same technique described in Sect. 14.11 We find that 293 of the 
24Ss are also detected in Ha (within 376). An inspection of the 
Ha image shows that most bright 24Ss with no Ha counterpart 
actually lie near an Hn region, with an offset which, though 
larger than the photometric errors, could still be accommodated 
in case of asymmetries of the Ha shell. This is the case, for 
example, of some sources along the spiral arms. Some faint 
24Ss sources might instead really be absent in the Ha map. 
Indeed, the percentage of 24Ss undetected in Ha gets higher 
close to the completeness limit of our survey. We will analyse 
this finding in more detail in a forthcoming paper. We note 
however that the luminosity function of the 24 Ss with Ha 
counterpart is totally compatible with the one for the whole 
sample presented in Sect. 16.31 This is because contamination of 
the 24 pm sample by non star forming regions happens mostly 
below the completeness limit of the derived luminosity function. 

To estimate the SFR in discrete sources, we adopted the same 
conversion methods used for the global SFRs and, in order to 
estimate the TIR luminosity following Calzetti's formula, we re- 
strict our sample of 293 sources to the 186 sources which are also 
detected in the IRAC bands. The resulting SFRs are compared 
in Fig. [T71 where the SFR is derived from Ha without any aver- 
age or individual extinction corrections and the FIR fluxes have 
been converted to SFR as in the previous subsection: SFR [M© 
yr" 1 ] = 8 x 1CT 44 L(FIR). The scatter about the line of slope 
unity is large, roughly a factor of 5 (l<x), and it is difficult to as- 
sess the presence of any systematic deviation. The scatter itself, 
given its magnitude, cannot be imputed to extinction only. Most 
of the variance is likely linked to the uncertainties intrinsic to the 
methods used to infer SFRs. A very similar plot can be obtained 
using the SFR estimated only from the 24 pm emission (see pre- 
vious subsection) of the discrete sources and is consequently not 
shown here. 

We can compare the IR colours of the latter sample with 
those of the complementary sample made up by the 121 24Ss 
with no Ha counterpart, selecting sources detected in all IRAC 
and MIPS-24 bands. One such diagram is in Fig. [18] and should 
be compared with the equivalent diagnostic diagram of Sect. [5] 
The 24Ss without Ha actually seem to be less dusty sources than 



those emitting also in Ha, a trend confirming that 24Ss with no 
Ha are not more embedded into dust. 



8. Summary and conclusions 

This is the first in a series of papers aiming to investigate the star 
formation history in M 33 by means of the IR data from Spitzer. 
We retrieved and reduced IRAC and MIPS data producing maps 
roughly the size of the entire star forming galaxy disk. These 
were used to study the nature of the MIR emission, both globally 
(large scale structure of dust emission) and locally (by means of 
discrete sources). Our main results are the following: 

- On large scales, the 3.6/4.5 and 5.8/8.0 ratios are remarkably 
constant. The 3.6 and 4.5 pm emission mainly arises from 
point sources while the 5.8 and 8.0 pm bands present more 
diffuse emission following the morphological pattern of the 
galaxy (spiral arms). The 8.0/24 ratio is also rather constant 
suggesting a tight link between PAH and VSG properties. 
On the other hand, the 3.6/5.8 and 4.5/8.0 ratios reveal the 
flocculent spiral structure in M 33, because of the increasing 
dominance of the ISM toward longer wavelengths. The dis- 
tribution of cooler dust is investigated via the MIPS images: 
the 24 and 70 pm emissions match closely, while the 160 pm 
band reveals a more diffuse character in the cold-dust emis- 
sion. The 24 and 70 pm are also more tightly correlated with 
the Ha emission. 

- In the four IRAC and 24 pm MIPS bands, we performed pho- 
tometry of catalogued discrete sources of various type, in- 
cluding Hn regions, SNRs and PNe. We present the detection 
rates in each band for the various types of object, together 
with their IR SEDs. For all types of objects, a drop at 4.5 pm 
is observed, likely due to the contribution of emission fea- 
tures or to the decreasing contribution of photospheric emis- 
sion from the central star. 

- The IR colours of Hn regions, SNRs and PNe are used to 
define diagnostic diagrams. We observe a continuity among 
the various types: the Hn regions appear more dusty than the 
PNe, which behave more like stellar objects. We also com- 
pared these colours with colours predicted by models (such 
as Cloudy, DUSTY, SB 99, pure blackbody emission). The 
diagnostics based on IRAC alone cannot disentangle the dif- 
ferent theoretical models, but usage of hybrid IRAC +24 pm 
colours allows the removal of some of the degeneracies. Hn 
regions are better modelled by pure Cloudy Hn regions (in- 
cluding PAHs) and by DUSTY models with inner-shell tem- 
peratures around 600 K and low extinction. 

- We have compared the colours of some of the known-type 
objects in M 33 with similar objects in our Galaxy (Hn re- 
gions) or in the LMC (SNRs) and, although with some scat- 
ter, we found good agreement for the location of the sources 
in the colour diagrams. 

- From our 24 pm image, we extracted a blind catalogue of 
515 discrete sources, on which we performed photometry in 
the IRAC, and 24 pm MIPS bands and in Ha. Their colours 
match closely those of known-type sources, and in particular 
those of Hn regions. Consequently, the models which better 
describe the nature of the 24 pm sources are also pure Cloudy 
Hn regions and low-extinction warm (Ti n = 600 K) DUSTY 
models. 

- The total 24 pm flux of the 515 discrete sources amounts 
to 1.4 x 10 41 erg s" 1 while the diffuse emission, that is the 
residual flux, amounts to twice this value. 
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- We present the IR luminosity function of our 24 jim sources, 
as a function of the 24 fim flux and as well as a function of 
the total IR luminosity (estimated from the 8.0 and 24 //m 
bands). The luminosity ranges from values typical of the 
faintest ionising stars to those of very large star formation 
complexes. The shape (double- sloped) is the same as that 
observed for Galactic Hn regions and star complexes, though 
AGB stars could contribute significantly to the faintest bins. 

- Within a radius of 5 kpc in M 33, we estimate a total IR lu- 
minosity of 1.0 x 10 9 L and the FIR luminosity is about 
5.8 x 10 8 L . This translates into an estimate of the SFR of 
0.2 M yr" 1 , in very good agreement with the SFR inferred 
from Ha emission. 

- The relation between 24 /mi emission and Hn regions is also 
investigated by searching directly for Ha counterparts of the 
24 jim sources. About half of the 24 jim sources are very 
closely associated with Ha while the remaining ones are 
more distant from an Ha knot than the allowance by posi- 
tional uncertainties. The SFR for the single sources in IR 
and Ha are consistent but the scatter is rather large and not 
imputable to extinction alone. 

- The IR colours of the 24 jim sources with and without Ha 
counterparts are quite similar. If anything, upon comparison 
with the theoretical diagnostics, those with Ha appear more 
dusty than the others. 
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Fig. 8. The position of catalogued sources with photometry in all 
four IRAC bands superimposed on the 4.5 //m IRAC image. The 
upper panel shows the IR selected (open squares), radio selected 
(pluses) and optically selected (open circles) Hn regions. The 
lower panel shows the locations of the PNe (open circles), op- 
tically selected (crosses) and radio selected (pluses) SNRs. The 
images are about 1° on a side, north is up and east to the left . 
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Fig. 9. Mean IR SEDs of the different kind of sources. Open 
symbols represent the Hn regions selected in IR (squares), radio 
(triangles) and optical (circles). The (blue) crosses represent the 
SNRs optically selected and the (blue) pluses the SNRs radio se- 
lected. The filled (red) diamonds represent the PNe. The (green) 
tripods represent the 24 /mi unknown sources. The grey vertical 
stripes depict the widths of the IRAC and MIPS-24 filters. 
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Fig. 10. Spectra obtained with Cloudy code, for the four mod- 
els. The parameter pairs in size-density space with density given 
first in cm -3 followed by outer radius in pc are represented by 
the (blue) short-dashed line for (100,10); the (cyan) dotted line 
(20,30); the (green) long-dashed line (10,60); and the (red) solid 
line (1000,2). 
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Fig. 11. The resulting DUSTY spectra. We plot only 5 values 
of A y , namely 0.1, 2.15, 10.0, 21.5, and 100 (shown as long- 
dashed, short-dashed, dotted, dot-dashed, and solid lines, respec- 
tively). Also shown (as a heavy solid line which ends at X ~ 
lOOyum) is the input stellar spectrum of a 3Myr starburst from 
SB99. The upper panels show Ti n = 300 K, and the lower ones 
T in = 600 K. The left panels show thick shells (R out = 10000), and 
in the right panels thin ones (R ou t= 100). 
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Fig. 12. Bottom panel: Blackbodies are depicted by filled 
(black) circles. NGC 7023 is represented by a (yellow) triangle, 
the ISM by Draine by an inverted (cyan) triangle, the SB 99 simu- 
lations by (grey) diamonds. The Cloudy Hn region is represented 
by filled (green) squares. Finally, the DUSTY models are repre- 
sented by filled (Ti n = 300 K) and empty (Ti n = 600 K) circles, 
for a variety of shell thickness Y ou t = 10 (solid), 100 (dashed), 
1000 (dash-dot) and 10000 (dotted line). Different values of A y 
(logarithmically from 0.1 to 100 mag) for a given Ti n and Y out 
are connected by the different lines; the highest value of Ay is 
marked by x. Middle panel: The open symbols represent the 
Hn regions selected in IR (squares), radio (triangles) and op- 
tical (circles). The (blue) crosses represent the SNRs optically 
selected and the (blue) plus the SNRs radio selected. The filled 
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Fig. 13. Same as Fig. El but for different IRAC flux ratios. 
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Fig. 14. Same as Fig.El but now including the MIPS 24 jum flux 
ratio. 
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Fig. 15. Diagnosti c diagrams in magnitude, for comparison with 
ICohen et aiT(l2006l) . The open symbols represent the Hn regions 
selected in IR (squares), radio (triangles) and optical (circles). 
The (blue) crosses represent the SNRs optically selected and the 
(blue) plus the SNRs radio selected. The filled (red) diamonds 
represent the PNe. Two radio Hn regions are lying outside the 
limits of the displayed plot. 
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Fig. 16. Bottom panel: Cumulative distribution of the 24 //m 
sources. The 24 /mi fluxes are reported F v (24 jim) (bottom axis), 
as well as the TIR luminosity L y (TIR) (top axis). The linear best 
fits are depicted by a (red) solid line for the F y (24 //m) fluxes 
between 1 and 70 mJy, and by a (green) dashed line for the 
F y (24 jim) fluxes greater than 70 mJy. Top panel: Cumulative 
distribution of the 24 jim sources detected also at 8.0 /mi, the 
TIR luminosity is calculated using for each source the effective 
8 yum/24 jim flux ratio. 
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Fig. 17. Comparison of the SFR calculated from the Ha lumi- 
nosity and from the IR emission based on the 8.0 and 24 fim 
emissions. The dashed line shows the SFR(Ha) = SFR(IR) rela- 
tion. 
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Fig. 18. Diagnostic diagram involving IRAC colours for the 
24 yum sources having a Ha counterpart (open circles) and for 
the 24 fim sources with no Ha counterpart (crosses). 



